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ABSTRACT 

We measure the rest-frame B-band luminosity function of red-sequence galaxies (RSLF) of 
five intermediate-redshift (0.5 < z < 0.9), high-mass (a > 950 km s _1 ) clusters. Cluster galaxies 
are identified through photometric redshifts based on imaging in seven bands (five broad, and 
two narrow) using the WIYN 3.5m telescope. The luminosity functions are well-fit down to 
Mg + 3 for all of the clusters out to a radius of i?200- For comparison, the luminosity functions 
for a sample of 59 low redshift clusters selected from the SDSS are measured as well. There is a 
brightening trend (Mg increases by 0.7 mags by z=0.75) with redshift comparable to what is seen 
in the field for similarly defined galaxies, although there is a hint that the cluster red-sequcncc 
brightening is more rapid in the past (z > 0.5), and relatively shallow at more recent times. 
Contrary to other claims, we find little evidence for evolution of the faint end slope. Previous 
indications of evolution may be due to limitations in measurement technique, bias in the sample 
selection, and cluster to cluster variation. As seen in both the low and high redshift sample, 
a significant amount of variation in luminosity functions parameters a and M* exists between 
individual clusters. 



Subject headings: galaxies:cluster: general — galaxies: evolution — galaxies: luminosity function 



1. Introduction 

In the nearby Universe, galaxy clusters are 
dominated by early type galaxies (Dressier 1980) 
along the so-called "red sequence," the name for 
which derives from the tight color-magnitude re- 
lation observed for these galaxies (Visvanathan & 
Sandage 1977). The most luminous cluster galax- 
ies, located at the tip of the red sequence, appear 
to have relatively little star formation since z ~ 3 
(Bower et al. 1992; Ellis et al. 1997; Stanford et 
al. 1998; Kelson et al. 2001). Observations of 
high redshift (z > 1) clusters are consistent with 
this inference, indicating at least the bright end 
of the color-magnitude relationship has only pas- 
sively evolved since z ~ 1.5 (Stanford et al. 1997, 
Mullis et al. 2005, Stanford et al. 2005). How- 
ever, evidence for recent (< 5 Gyr) bursts of star 
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formation in today's lower-mass cluster galaxies 
(Poggianti et al. 2001, Conselice et al. 2003) 
indicates not all cluster galaxies have the same 
homogeneous star formation histories. While the 
luminous red sequence in clusters may have long 
been in place, the extension of the red sequence to 
lower luminosity may be a relatively recent phe- 
nomenon. 

Recent, deep imaging programs have led to the 
claim of a deficit of red sequence at low luminosi- 
ties in high redshift clusters (Nakata et al. 2001, 
de Lucia et al. 2004, Goto et al. 2005, Tanaka et 
al. 2005). This would suggest many of the galaxies 
seen along today's red sequence may be late addi- 
tions, perhaps culled from a quenched and faded 
blue population. This fits neatly with a notion 
that the red sequence is not monolithic, as evinced 
by the seminal studies of the Virgo cluster lumi- 
nosity function (Sandage, Binggeli & Tammann 
1985). However, the deficit may not be ubiqui- 
tous (e.g., Andreon 2006). Differences in the anal- 
ysis methods, clusters samples, and the inherent 
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limitations in different data sets (i.e., the depth 
of available data) have left the observational evi- 
dence for a deficit inconclusive. 

Differences between studies are also compli- 
cated by the possible existence of evolutionary 
phenomenon that is a function of cluster mass. 
Substantial differences exist between the cluster 
and field red sequences (de Propris et al. 2003, 
Croton et al. 2005), and there is evidence for 
differences in the luminosity function with cluster 
mass, both locally (Hansen et al. 2005, Hilton et 
al. 2005) and at intermediate redshift (Koyama et 
al. 2007, Gilbank et al. 2007). However, evidence 
to the contrary has also appeared in the literature 
(Barkhouse et al. 2007, Andreon 2007). Given 
these disparate claims, it is difficult to assemble 
a clear picture of the red-sequence evolution - a 
challenge again exacerbated by differences in clus- 
ter samples and measurement techniques between 
studies. 

A common theme repeated in many of the stud- 
ies is the large cluster to cluster variation seen in 
any of the parameterizations of the cluster popu- 
lations, particularly the luminosity function. Es- 
timates of the dwarf-to-giant ratio (DGR, de Lu- 
cia et al. 2006) and Schechter-function (Schechter 
1976) fits to the luminosity function (Barkhouse 
et al. 2007, Andreon 2007) both show significant 
cluster-to-cluster variation. Unique features also 
have long been identified in the shape of lumi- 
nosity functions of local clusters (Biviano et al. 
1995, Yagi et al. 2002), indicative of a multi- 
component population even among the red-galaxy 
population. Traditionally, deep studies of the so- 
called red-sequence luminosity function (RSLF) in 
local clusters have fit two functions (cither double 
Schechter functions or a Gaussian and Schechter 
function) to the distribution of galaxies (Sandage 
et al. 1985, Jerjen & Tammann 1997, Popesso 
et al. 2005). The bright end of the luminosity 
function is well described by a functional form 
with a sharp turnover at lower luminosities, while 
the fainter luminosity-function component is gen- 
erally found to be rapidly rising at lower lumi- 
nosities. Variations in the relative amplitude of 
these components may well drive overall variations 
in observed cluster luminosity functions. Studies 
of intermediate redshift clusters are dominated by 
this brighter component, due to depth limitations. 
While such studies rarely reach deep enough to 



well-characterize the rapidly-rising faint end, its 
modulation may significantly impact the observed 
count-parametrization. 

In this paper, we explore the red-sequence 
population in five high-density regions contain- 
ing massive, bona fide, and well-studied clusters 
between 0.5 < z < 0.9. Our study is based on 
deep, multi-band imaging data from the WIYN 
3.5m telescope and extant spectroscopic data. The 
depth of our WIYN data is comparable to other 
cluster studies at similar redshifts, and allows an 
independent assessment of the cluster RSLF sig- 
nificantly below the knee in the bright end of the 
luminosity function. The extent to which this, or 
any other extant study at intermediate redshift, 
can address the faint-end of the cluster RSLF 
is a point we address. For comparison, we in- 
clude measurements of the luminosity function for 
a large sample of low redshift clusters from the 
Sloan Digital Sky Survey measured in the same 
manner as our intermediate redshift sample. 

In §2 we summarize the observation data and 
highlight our measurement techniques. In $3l we 
present the luminosity functions measured for each 
of the clusters as well as the tests to confirm and 
assess the reliability of the measured luminosity 
functions. In £|3 . 31 the low-redshift cluster sample 
is described. Finally, we discuss the observed evo- 
lution in our luminosity functions (§4), compare 
to results in the literature, and discuss the impli- 
cations of our results in §5. Throughout this work, 
we adopt H = 70 km s _1 Mpc -1 , flMatter = 0.3, 
and f^A = 0.7; all absolute magnitudes are in the 
rest-frame (Johnson) B-band in the Vega system. 

2. Observations and Analysis 

Observations were obtained between 1999 Oc- 
tober and 2004 June, with the WIYfsQ 3.5m tele- 
scope's Mini-Mosaic Camera (0.14" per pixel and 
9.6' x 9.6' field of view, 0.85" FWHM median see- 
ing) as part of the survey described by Crawford 
et al. (2006) and Crawford (2006). Cluster images 
were taken in the Harris UBRI, Gunn z (Schnei- 
der, Gunn, & Hoessel 1983), and two narrow band 
filters. The narrow band filters (A/AA ~ 70) were 
specifically designed to detect [OIIJA3727 and ad- 
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jacent continuum at each cluster redshift. Table 
Q] lists the five clusters observed, and their salient 
attributes. 

Data analysis proceeded on deep mosaic images 
created from reduced Mini-Mo frames, flat to 1% 
of their initial sky values. The deep mosaics were 
photometrically calibrated through a variety of in- 
dependent methods yielding uncertainties below 
2%. Object detection was performed on the ri- 
band (typically deepest) images using SExtractor 
(Bertin & Arnouts 1996) with the criterion that 
objects contained > 20 contiguous pixels above 
3<r of the sky noise. This corresponds to a min- 
imum signal-to-noise (S/N) of 13.4 in an equiva- 
lent circular aperture of 0.35 arcsec radius (just 
slightly smaller than the typical seeing disk) . The 
S/N at the 50% detection limit (see below) will 
be higher since most objects intrinsically at the 
specified limiting size and surface-brightness will 
have some fraction of the pixellated signal noise- 
aberrated below the detection threshold. As the 
S/N depends on the aperture size and tends to 
peak around the the half-light radius, it is essen- 
tial to define the S /N used in any analysis to avoid 
confusion. 

Two types of magnitudes were used for all cal- 
culations. Seeing-matcheoH aperture magnitudes 
with a radius corresponding to 7.5 kpc at each 
cluster redshift were used for all colors and pho- 
tometric redshifts. This is equivalent to aper- 
ture diameters between 1.93 and 2.45 arcsec, or 
roughly 2.5x stellar FWHM on average. This 
aperture defines the S/N value for the R-band that 
we use in subsequent analysis, particularly in the 
next two sub-sections concerning detection com- 
pleteness and photometric-redshift errors. Tai- 
lored magnitudes (Crawford et al. 2006) based 
on the curve of growth, concentration index, and 
Petrosian ratio are used for total apparent and 
absolute magnitudes. Absolute magnitudes and 
colors were calculated for all objects using photo- 
metric (see below) or spectroscopic redshifts. K- 
correction calculations adopt method 4 of Ber- 
shady (1995). 
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2.1. Detection Completeness 

We determine the completeness limit of our 
data through Monte Carlo simulations similar to 
what is described in Bershady, Lowenthal & Koo 
(1995). In each detection image, we extract a sam- 
ple of bright, representative objects of different 
sizes. Objects are dimmed, and then re-inserted 
into the image. For each half-magnitude interval, 
we insert 50 objects randomly across the usable 
field of view, repeat the detection procedure, and 
repeat the process 10 times, and for a minimum of 
three different object sizes. For small, unresolved 
objects, the 50% completeness magnitude is typi- 
cally around R ~ 25.5 magnitudes. For objects 
with lower surface brightness, our completeness 
limits can be significantly (~ 1 mag) more shallow. 
The R-band completeness curves are presented in 
Figure [T] For guidance as to the depth of these 
data, we mark in the Figure the location of M* 
at the appropriate cluster redshift based on the 
evolving RSLF estimated for the field (Willmer et 
al. 2006). When correcting for incompleteness, we 
calculate the correction relative to the measured 
size of each object. S/N in seeing- matched aper- 
tures at the 50% and 90% detection completeness 
are roughly 20 and 25, respectively, which is well 
above the limiting S/N for detection. 

For luminous red galaxies, the completeness 
limit is relatively shallow due to the diffuse profile 
shape and large size of such galaxies. Although a 
significant portion of the light is concentrated into 
the central regions, most of the light is in the low 
surface brightness wings of the galaxy. However, 
these systems are not at the detection-limit of our 
survey at the cluster redshifts. Instead, intrinsi- 
cally faint (low luminosity), red galaxies are sam- 
pled near our detection limits, for which the above 
detrimental properties are not in play. Three nat- 
ural phenomenon are working in our favor in terms 
of detecting these lower-luminosity galaxies. First, 
the light profiles of lower-luminosity red-sequence 
galaxies are closer to an exponential distribution 
than an r 1 / 4 (Graham & Guzman 2003), thereby 
enhancing detectability. Second, according to the 
luminosity-size distribution of red galaxies (Shen 
et al. 2003, Mcintosh et al. 2005), we expect the 
faint, red galaxies to be unresolved in our data set, 
and hence our deepest detection limits are rele- 
vant. We return to, and illustrate these two points 
later in §3. Finally, lower-luminosity galaxies are 
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expected to be slightly bluer than their luminous 
red-sequence counterparts due to a slope in the 
color-magnitude relationship. Hence, the effect of 
the k-correction for these galaxies at the R-band 
detection limit is less severe. On this basis, we 
compute the depth of our survey relative to the 
50% completeness limit for unresolved objects in 
that field. These depths are listed in Table HJ and 
relative to the field RSLF, they range between 2.6 
and 3.6 mag below M*. 

2.2. Photometric Redshift Method, Accu- 
racy and Precision 

We have derived photometric redshifts for all 
objects with the highest precision possible. The 
photometric redshifts were determined through a 
hybrid of the template and training-set methods 
similar to Csabai et al. (2003). A template grid 
of model spectral energy distributions (Bruzual 
& Chariot 2003), and their corresponding broad- 
and narrow-band fluxes, was created. The spectra 
cover a range of star formation rates, ages, and 
redshifts. Using objects with spectroscopic red- 
shifts and our observed multi-band fluxes, the grid 
was constrained (more specifically, morphed or ad- 
justed in flux space) to match the measured data. 
Over 500 objects with existing spectroscopic red- 
shifts in our combined survey fields (Dressier et al. 
1997, Ellingson et al. 1998, Postman et al. 1998, 
van Dokkum et al. 1999, Moran et al. 2007, Tran 
et al. 2007) trained our set of templates in this 
way. The method was tested with the following 
procedure. A single spectroscopic object was re- 
moved from the catalog. The grid was trained us- 
ing the remaining objects and then applied to the 
removed objects. The procedure was repeated so 
that an independent measure of the photometric 
redshift (and its error) was made for every spec- 
troscopic object. In Figure [2] and [3l we compare 
the spectroscopic and photometric redshifts mea- 
sured via this procedure. High quality photomet- 
ric redshifts can only be obtained with large val- 
ues of signal to noise as shown in Figure [H Pho- 
tometric redshifts are calculated with a precision 
cr 2 = 0.03 and an accuracy of S z < 0.01 for red ob- 
jects with S/N > 30 in the R-band. Blue objects 
(defined below in §2.3) performed slightly worse 
with cr 2 = 0.05 and 5 Z = 0.01 with S/N > 30. 
Both types of objects have much larger scatter for 
galaxies with S/N < 30. At the detection thresh- 



old, the typical scatter for red objects is expected 
to be approximately a z — 0.05. 

2.3. Impact on Luminosity Function Cal- 
culations 

Extensive simulations of the effect of photo- 
metric errors on the distribution of measured 
photometric-redshift errors, as well as the sub- 
sequent effect on the cluster-galaxy luminosity 
function were performed. The latter comes about 
because we define cluster membership based on 
photometric redshift (§3) when spectroscopic red- 
shifts arc not available. Photometric redshift er- 
rors at high signal to noise (S/N > 10) are small, 
and have a Gaussian distribution. At low signal 
to noise, however, the error distribution is a com- 
bination of a Gaussian core around the fiducial 
redshift plus a catastrophic error component due 
to near-degeneracies in multi-color space. Com- 
bined with our definition of cluster membership, 
the photometric-redshift errors cause a significant 
portion (> 10%) of the cluster population above 
our detection limit to be identified as field galax- 
ies beyond our redshift selection window. Because 
of large cluster over-densities, contamination from 
the field into the cluster sample is minimal. How- 
ever, the opposite is true when we construct the 
field-population estimate for the RSLF, as we 
show below in §2.3. Consequently, the nature 
of the correction is to account for a net loss, or 
incompleteness in the bona-fide number of cluster 
galaxies, and vice-versa for the field population. 
These corrections are independent of our detection 
completeness except in so far as both are corre- 
lated with apparent magnitude (or S/N). Statis- 
tical corrections, based on our simulations, are 
applied to our luminosity functions as a function 
of apparent magnitude. At the completeness limit 
in the images, the average correction is 20% with a 
maximum of 40% for CL0016. For the field-galaxy 
luminosity function that we construct below as a 
quality check, we push these corrections farther to 
illustrate the quality of these corrections. In §3.2.1 
we discuss tests of the above issues of complete- 
ness and contamination in the specific context of 
our cluster selection function. 
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2.4. Quality Check: The Field RSLF at In- 
termediate Redshifts 

To demonstrate our methodology for selecting 
and computing the statistical distribution of red 
cluster galaxies, we first construct the field RSLF 
luminosity function, derived from our data, in 
two redshift bins that should be uncontaminated 
by rich, over-densities, i.e., foreground and back- 
ground regions in the cluster fields (Figure[5]). The 
low redshift bin (z — 0.3 — 0.55) contains the 
imaging from fields MS1054, C11322, and C11604, 
whereas the high redshift bin (z = 0.7—0.95) is cal- 
culated based on data from C10016 and MS0451. 
In addition, one empty "field" image centered at 
13:24:50. lh +30:11:18.5°, which contains no evi- 
dence of any massive over-density, was included in 
both bins. Red galaxies were identified following 
the same prescription as Willmer et al. (2006): 

U -B < -0.032 * (M B + 21.52) + 0.204. (1) 

This definition, used throughout the paper, is 
based on a fit to the color-magnitude relation plus 
a -0.25 mag shift in zeropoint. 

Corrections for detection incompleteness and 
contamination due to photometric redshift errors 
(as described in the previous section) have been 
applied. The data are only presented to a mag- 
nitude limit where the completeness is > 50%, 
but in some cases the contamination corrections 
are larger. At low S/N (faint magnitudes), cluster 
sources with large photometric-redshift errors can 
contribute substantially to the field counts. Akin 
to the Eddington effect (Eddington 1913) where 
fainter, more plentiful objects have a more sig- 
nificant influence on brighter bins with less ob- 
jects, the majority of field counts at faint mag- 
nitudes can be due to contamination from clus- 
ter over-densities. For this reason, and to keep 
contamination corrections manageable, we calcu- 
late the luminosity function after excluding red- 
shift bins where, for a given luminosity, the con- 
tribution from the cluster is expected to be greater 
than 5%. For example, at R = 25 (Mb = -18.6) 
in the low-z sample, the error on the photomet- 
ric redshift is expected to scatter up to 10% of 
the galaxies from MS1054 into the redshift bin of 
Az = 0.5 — 0.6. Hence, in this magnitude bin, we 
would exclude galaxies with photometric redshifts 
between z = 0.5 - 0.6 in the MS1054 field from 
the low-z luminosity function while still including 



galaxies with lower redshifts. An appropriate vol- 
ume would be calculated for this magnitude bin 
according the corresponding redshift window. 

In Figure O we compare our raw and corrected 
field RSLF to DEEP2's Schechter-function param- 
eterization of the same (Willmer et al. 2006). 
A reduced-% 2 statistic of this parameterization 
with respect to our data-set indicates a reason- 
able match between results, despite our smaller 
survey volume and use of photometric redshifts. 
We reach a similar conclusion comparing our field 
blue-galaxy luminosity function with correspond- 
ing results Willmer et al. (2006). This is inter- 
esting because in the case of the blue-galaxy lumi- 
nosity, the field sample we construct has smaller 
contamination errors. On balance, these results 
indicate we are able to control our corrections 
for incompleteness and photometric errors reliably 
down to and below the level we will use in our clus- 
ter luminosity function analysis. 

3. Cluster Red-Sequence Luminosity Func- 
tions 

3.1. WIYN Intermediate Redshift Sample 

Red-sequence cluster galaxies were identified 
with (i) the same color-luminosity relation (Eq. 1) 
as for the field, and (ii) as having redshifts within 
±0.05 of the cluster redshift. Both of these se- 
lection criteria for our data can be seen in Fig- 
ure [5] Most of the clusters do exhibit a strong 
red sequence even when data are included from 
the entire field of view, which is much larger than 
the cluster core-radius. All fields show an over- 
density at the redshift of the cluster. For those 
two clusters (C11322 at z=0.75 and C11604 at 
z=0.9) which do not show a strong red excess in 
the color-redshift plot due to contrast (Figure 6, 
right panel), they still show a well-populated red- 
sequence in the color-magnitude diagram (Figure 
6, left panel). The raw histogram in absolute mag- 
nitude of red-sequence galaxies, normalized by the 
selection volume, are displayed in Figure [7] (left 
panel) for a selection radius of i?20(H- (The volume 
was estimated for each cluster from the selection 
in redshift and radius.) Since the cluster is ex- 
pected to cover a much narrower range in redshift, 
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the volume, while well-defined, is an over-estimate. 
Aside from this volume normalization, this his- 
togram represents the uncorrected measurement 
of the red-sequence luminosity function in terms 
of shape and luminosity normalization. 

Several corrections were applied to the data to 
account for observational deficiencies: 

(a) The first correction was for detection com- 
pleteness, computed from the Monte Carlo simula- 
tions (§2.1). The correction is calculated for each 
galaxy based on its apparent magnitude and size. 

(b) Because we sample a much larger volume 
than the cluster, a portion of the counts are ex- 
pected to be from the field around the cluster. We 
therefore subtract an appropriately volume-scaled 
field sample based on the RSLF of Willmer et al. 
(2006), which we note is always a small fraction of 
the total uncorrected cluster counts (e.g., < 10% 
within 1 Mpc). 

(c) Finally, the incompleteness due to photo- 
metric redshift errors was applied (§2.2), again 
based on simulations. Specifically, we computed 
the fraction of galaxies as a function of apparent 
magnitude that will be missed (scattered out of 
our redshift selection window). We can ignore as 
negligible the fraction scattered into our redshift 
selection window from the field, because to first 
order this cancels with the number of field galax- 
ies in our volume scattered out, and further we 
subtract a field component from our volume, as 
noted above in (b). By assuming a single spectral 
energy distribution for the red galaxies, this cor- 
rection can then be applied to the binned counts 
as a function of absolute magnitude. 

Absolute magnitude bins containing galaxies 
beyond our specified completeness limits were 
excluded from our calculations, and our pri- 
mary results are quoted down to 50% detection- 
completeness limits. The final, corrected data 
within i?2oo for these limits are presented as solid 
points in Figure [7] (left-hand panel) . We also use 
90% detection-completeness limits in some cases 
to further demonstrate the robust nature of our 
statistical results. These are indicated by the ver- 
tical dotted lines in this Figure. 

For each cluster, we calculate the best-fit 
Schechter-function to the luminosity distribution 
by varying a, M* , and $* to find the smallest 
reduced x 2 . Because of the nature of our vol- 



ume normalization, we focus here on a and M* , 
not $*, with our primary scientific emphasis be- 
ing on a. We plot (Figure El right) the error 
ellipses calculated from x 2 measurements based 
on the best fit model. Again, these are for data 
within -R200 projected about the cluster center. 
The luminosity function also was calculated over 
several different selection radii: 1 Mpc, 0.5i?2oo, 
and 0.25i?2oo- These results and those for i? 2 oo 
are shown in Figure [5] and listed in Table O 

3.1.1. Trends with selection radius 

Substantial differences in the shape of the lumi- 
nosity function are seen at smaller selection radii, 
consistent with local clusters (Lobo et al. 1997, 
Popesso et al. 2006), and presumably due to a 
morphology-density relation in the dwarf-giant ra- 
tio. This effect is illustrated in our study in Figure 
8. Our clusters show a general trend of a flatter 
(a ~ —1) luminosity function with increasing clus- 
ter radius, albeit with significant scattered espe- 
cially for C11322 and C11604, where our errors are 
largest. In the literature, Lobo et al. (1997) find 
a steeper faint-end slope in the central regions of 
Coma as compared to groups around the outskirts. 
Popesso et al. (2006) found that the brightness of 
the faint-end luminosity function increased with 
increasing cluster radius in the context of dou- 
ble Schechter-function fits. For a single Schechter- 
function fit, this would result in a flatter [a ~ —1) 
fit to the luminosity function for larger selection 
radius, which is seen in our cluster sample. RSLF 
shapes are far more uniform with a i?2oo selection, 
but field contamination becomes much greater es- 
pecially for lower-density clusters. For direct com- 
parison with other intermediate redshift work this 
result indicates it is critical to make comparisons 
within the same selection radius, preferably rela- 
tive to -R200 for each cluster. This is a point we 
return to later. In general, however, we are obliged 
to use the R = 1 Mpc aperture. For our clusters, 
this is between 0.40-0.66 i?200- 

3.1.2. Comparison to literature 

A comparison can be made between our luminosity- 
function measurements and four previous mea- 
surements for three clusters. 

C10016. Tanaka et al. (2005) find M* B = -20.14± 
0.6 (assuming (B - V) = 0.9) and a = -0.64 ± 
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0.4. for C10016 in a selection radius that is close 
to -ffeoo- Relative to our measurement, this is 
+1.2 ± 0.7 mag fainter in M B (a modestly sig- 
nificant difference), but only +0.2 ± 0.4 shallower 
in a. In actuality, however, Tanaka et al. use a 
local density criteria for their selection. Given the 
well-known elongation of the central, high-density 
region of C10016 (e.g., their Figure 6), this results 
in an effective selection radius which is substan- 
tially (~ 50%) smaller. Indeed, comparing to our 
0.5i?2oo selection yields very similar values (see 
Figure 8), with their M* B only +0.2 ± 0.7 mag 
fainter, and their a only +0.042 ± 0.4 steeper. 
These are well within the lcr measurement errors, 
and small in an absolute sense as well. 
MS 1054- Our measurements for MS 1054 agree 
closely with those of Andreon (2006) for his 1 
Mpc selection radius: M* B = -21.0 ± 0.2 and 
a = -0.80 ± 0.12, which is 0.3 ± 0.4 different in 
M* B and +0.04 ± 0.24 different in a. 

However, Goto et al. (2005) find a value of 
M* B = -21.15 ± 0.22^ and a = -0.09 ± 0.27 for 
objects with ^775 — zgso > 0.5 within a selection 
radius of R — 1.35 Mpc from a spectroscopic sam- 
ple. Our value for M B — —21.16 within -R200 
(1.82 Mpc) is negligibly different than theirs, but 
we find a = —0.58. As Andreon (2006) points 
out, the spectroscopic completeness in Goto et al. 
at the faint end is only 20%, which may lead to 
a severe turnover. Furthermore, they use a sin- 
gle color cut to identify red galaxies instead of 
a luminosity-dependent color-cut, as done in this 
work. Consequently, bright red sequence galaxies 
are likely to be included whereas faint red galax- 
ies are likely to be preferentially excluded due to 
signal-to-noise considerations and the slope of the 
red sequence. If not properly accounted, this effect 
could be manifest as an apparent 'deficit' of faint 
red galaxies compared to bright ones. Their color 
also does not span the Balmer Break at z — 0.83. 

As a further check of the Goto et al. mea- 
surements, we note the luminosity function they 
derive for early-type galaxies is inconsistent with 
their measurement for red galaxies. Early-type 
galaxies are selected in their study via a visual 
(qualitative) morphological classification based on 
HST images. On the other hand, their value of 



4 The value from Goto et al. (2005) is converted from 1775 
AB magnitudes 



a = —0.54 ± 0.13 for early-types galaxies is very 
comparable to our measurement. The early- type 
galaxies at this redshift are still overwhelmingly 
red (and dominate the counts of red galaxies in 
MS1054; van Dokkum et al. 1999). No population 
of bright, red, late-type galaxies exist in MS1054 
that would bias the luminosity function in such a 
manner as might be inferred by the Goto et al. 
measurements. Based on the luminosity-function 
agreement between studies, we conclude their mor- 
phological selection is more reliable, in this case, 
than their color selection. 

C11604-. Andreon (2008) measured the luminosity 
function within the central 0.45 Mpc of the cluster 
from two-band HST imaging. He derived a value 
of a = —0.67 + 0.33 for the slope of the faint end. 
Our value of a = —0.86 ± 0.46 within an aperture 
of R=0.37 Mpc agrees within the lcr limits. 

3.2. Reliability of the Intermediate-Redshift 
RSLF Determination 

To further verify the accuracy of our results, we 
performed a number of tests of our data, reported 
here, including four tests of the measurement re- 
liability given the known amplitude of random er- 
rors, two tests investigating the systematic impact 
of the selection function for red galaxies, and a fi- 
nal test of the surface brightness bias present in 
our sample. Each of these tests validate a differ- 
ent aspect of our cluster measurement and reveal 
the quality of our data and analysis, and the ro- 
bustness of our results on the RSLF measurement. 

3.2.1. Random- Error Effects 

The first test we conducted was to measure the 
completeness and reliability of detecting red galax- 
ies within the cluster volumes selected. We calcu- 
lated the colors of a red galaxy (from synthetic 
spectra) at the redshift of each cluster, and at a 
redshift ±0.1 and 0.15 about the cluster redshift, 
i.e., at 2 and 3 times the distance in redshift as 
our nominal selection cut (± 0.05 about the clus- 
ter redshift). We then simulated the appropriate 
multi-band magnitudes and errors for the object 
for R = 19 to 26, in steps of 0.5 mags. For every 
magnitude interval and redshift bin, the simulated 
galaxy's photometry was realized 100-500 times as 
a perturbation about the nominal colors, draw- 
ing statistically on the error distribution for each 
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band to determine the perturbation. For each re- 
alization the photometric redshift was calculated. 
To determine completeness we simply counted the 
fraction of the galaxies simulated at the cluster 
redshift that retained photometric redshifts within 
± 0.05 of the cluster redshift. To determine the 
reliability, we calculated the percentage of galaxies 
simulated at ± 0.1 and 0.15 away from the cluster 
redshift had photometric redshifts estimated to lie 
within ± 0.05 of the cluster redshift. Such galaxies 
would be selected in our scheme as red-sequence 
cluster galaxies. This percentage was then normal- 
ized by the ratio of the number of galaxies in the 
cluster vs. the number of galaxies in the volume 
from z c i uster ± 0.15, as estimated in the real data. 
The contribution from the more distant shell was 
vanishingly small compared to the nearer shell, in- 
dicating our simulation should be quite accurate. 

The results of this simulation are shown in Fig- 
ure 9. From this it is clear that contamination 
is always <10%, and typically only a few percent 
even at the 50% source-detection limit. The se- 
lection completeness is 0.78 ± 0.15% at the 50% 
detection limit, and 0.86±0.13% at the 90% detec- 
tion limit. This is roughly what one would expect 
given the S/N and associated photometric-redshift 
errors at these limits. If anything, the results are 
somewhat optimistic, which can be understood in 
terms of the idealized nature of the simulation, 
e.g., the simulated galaxies are drawn from the 
same set that are used to detect and derive photo- 
metric redshifts. Overall, however, the simulation 
shows we are able to recover close to the expected 
number of cluster sources in a controlled situation 
resembling the actual data under analysis. 

The remaining three tests directly probe the de- 
rived luminosity function itself. 

In the second test, we exclude galaxies from our 
RSLF calculation where the detection-completeness 
was less than 90% (instead of the 50%). While 
we expect using corrected data down to the 50% 
completeness-limit is reliable (due to the exten- 
sive completeness simulations performed and the 
corrections derived therefrom), the robustness of 
our results is most directly shown by examining 
the truncated data set. We follow the same steps 
to calculate the luminosity function parameters 
as described in §3 but with fewer points. As can 
be seen by comparing the open to filled diamonds 
in Figure [7] (right), none of the clusters show a 



significant change in the parameterization of the 
luminosity function. At 90% completeness, all of 
our data extends beyond M* + 2. 

The third test was to create 100 realizations of 
each cluster RSLF from the measured errors in 
the corrected counts of galaxies identified as red- 
sequence cluster members. For each cluster, the 
observed luminosity function was convolved with 
the errors at each magnitude. Then, the param- 
eterization of the luminosity function was mea- 
sured for each realization, and the averages for 
the Monte-Carlo simulation were computed. The 
values found for M* and a are plotted as gray- 
scale in the left-hand plots in Figure [7] For all the 
clusters, the averages are well within the la er- 
ror measurements for the parameterization of the 
cluster luminosity function and the variance is of 
the same order as well. 

The fourth and final test was a more extensive 
test of the entire process of measuring the cluster 
luminosity function. Following the procedure of 
Toft et al. (2004), we produced ten realizations 
of our photometric catalogs for each cluster. In 
these realizations, the measured aperture photom- 
etry was smeared, in a statistical fashion, drawing 
from the photometric error distribution for each 
flux measurement. This means that the error dis- 
tribution for each realization is roughly larger 
than the initial (measured) catalog, although we 
did not update the effective error distribution for 
these realizations. With the smeared photometry, 
we recalculated the photometric redshifts. These 
new redshifts were then applied to the calculation 
of rest-frame properties and the selection of clus- 
ter galaxies. The luminosity functions were built 
using the same corrections and procedures as pre- 
viously described. Finally the parameterization of 
the luminosity function was measured for all of 
the clusters. The results of each of the realiza- 
tions are presented in Figure [7] For all but one of 
the clusters, the realizations produce results which 
show no significant difference in the mean from 
the measurement of the luminosity function. The 
dispersion in the luminosity function parameters, 
however, is larger, as expected from the additional 
noise introduced in the simulation process. For 
one of the smallest cluster in our sample (C11322), 
the results indicate a flatter (a ~ —1) luminosity 
function then measured in the single data alone. 
However, the original measurement is contained 
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within the spread of slopes that are found, and is 
not statistically significant. 

3.2.2. Systematics with Color- Selection 

To test what impact the specific selection of red 
cluster galaxies has on the derived RSLF, we se- 
lected galaxies from the ten realizations of each 
cluster catalog by applying perturbations to the 
color, magnitude, redshift, and radial selection 
functions. 

Changes to the magnitude zero-point of the 
color-magnitude relation made no changes to the 
selection of red galaxies due to the steep nature 
of the relationship. Large changes (> 0.1 mag) to 
the color zero-point did result in a small shift in 
the faint end slope of the luminosity function, with 
a tendency to find steeper downturns (more posi- 
tive a by +0.15) with a redward shift in the cut- 
off, and flatter slopes (more negative a by —0.14) 
with a blueward shift. This is qualitatively con- 
sistent with the canonical picture that the red- 
galaxy population and blue-galaxy population are 
characterized, respectively, by shallow and steep 
faint-end slopes to their luminosity functions. A 
similar trend was observed with changing the slope 
of the color-magnitude relationship. However, nei- 
ther change to the color selection-function would 
result in a measurement outside of the 95% con- 
fidence limits for the value found for the original 
color selection. This result is in qualitative agree- 
ment with previous analysis (Andreon et al. 2006, 
de Lucia et al. 2007, Barkhouse et al. 2007). 

3.2.3. Systematics with Redshift Window and 
Radius 

Variation in the redshift window do result in 
a small systematic change in the luminosity func- 
tion. The magnitude of the change, however, is 
much smaller than the errors on the measurement 
of an individual cluster: Aa < ±0.1 for a factor of 
two change in the redshift window. When the red- 
shift window is increased, the luminosity function 
becomes steeper (closer to a flat faint-end with 
a ~ —1), and becomes shallower (closer to a ~ 0) 
when the window is decreased. 

The systematic change in the luminosity func- 
tion with the change in selection window may 
indicate an underestimation of the photometric- 
redshift error-correction or contamination from 



field sources from the larger volume being investi- 
gated. Without extensive spectroscopic redshifts 
at faint magnitudes, it is difficult to conclude the 
source of this bias. Alternatively, this may be sim- 
ilar to the general trend seen in the change of slope 
with selection radius, albeit washed out by errors 
in photometric redshift, i.e., a steepening of the 
slope as the core is more preferentially sampled. 
Regardless, the small magnitude of the change 
provides confidence in our measurement and the 
robustness of our corrections. 

Small changes to the radial selection function 
do not result in significant changes to the luminos- 
ity function. Trends with larger variations in the 
radial selection, already noted, will be discussed 
again later. 

3.2.4- Surface- Brigtness Selection Effects 

Finally, we check that the depth of our obser- 
vations - in terms of surface-brightness sensitiv- 
ity - are sufficient to detect cluster members at 
these redshifts. As can be seen in Figure [101 the 
size-magnitude locus of the red sequence sources 
for most of the clusters is well within our surface- 
brightness detection limits. Only for the highest 
redshift cluster, C11604, are we truly in danger of 
missing some of the objects. 

However, if we are missing a significant amount 
of objects at the faint end, this will have the effect 
of causing the slope to fall more steeply than it 
should (a too large), which would mimic the as- 
trophysical effect claimed by others to be an evo- 
lutionary phenomenon. We do not see this effect 
in C11604 or MS1054, which both have slopes on 
par with our other clusters, and in fact C11604 has 
nominally the most negative a in our sample. 

3.3. The RSLF at z ~ 0.1 

A number of studies have measured the z ~ 
luminosity function from a variety of different 
sources using a wide range of techniques. To pro- 
vide a single, simple comparison to our body of 
work, we have measured the RSLF for a large sam- 
ple of Abell clusters (Abell, Corwin, & Olowin 
1989) based on SDSS imaging data and spec- 
troscopy (Adelman-McCarthy et al. 2008) using 
a similar procedure as for our clusters. Clusters 
are selected from the Abell catalog that also ap- 
pear in the SDSS DR6 imaging and spectroscopic 
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catalogs. For all the clusters that do appear in 
the SDSS data, we measure the velocity disper- 
sions of the clusters based on the spectroscopic 
data. From these measurements, we only include 
clusters with a minimum of 30 spectroscopically 
confirmed members. Clusters with extremely high 
velocity dispersions (a v > 1500 km s _1 ) or with 
significant difference in richness and a v are exam- 
ined individually to confirm that they are not close 
superpositions of two smaller cluster. Superpo- 
sitions are eliminated from the sample. Finally, 
cluster velocity dispersions are compared to val- 
ues found in the literature from Struble & Rood 
(1991), Wu, Xue, & Fang (1999), Miller et al. 
(2005), or Popesso et al. (2005). Clusters with 
large disparities between their literature value and 
that measured by the SDSS data were removed as 
well. Our final sample of 59 clusters have redshifts 
between 0.035 < z < 0.144 and a v > 500 km s _1 
with average values of z — 0.078 and a v = 886 km 
s -1 . Twenty-six of the clusters have a v > 900 km 
s _1 . The selected clusters are listed in Table [3] 

Using the SDSS photometric data, we have 
identified and analyzed the data following the 
same procedure as our cluster sample. The pri- 
mary difference is that we only use the five SDSS 
bands that are available and no narrow band 
data. For magnitudes, we used corrected Pet- 
rosian magnitudes following the recipe in Graham 
et al. (2005), and for colors, we use the SDSS fiber 
magnitudes, which are analogous to our aperture 
magnitudes. We have employed the same photo- 
metric redshift technique with a training sample 
created from the SDSS spectroscopic sample, cal- 
culated absolute magnitudes and rest-frame col- 
ors, and selected red sequence galaxies using the 
same selection function except shifted to z=0.08 to 
account for luminosity evolution. Cluster galax- 
ies were selected by having photometric redshifts 
within 0.05 of the cluster redshift. The luminosity 
function was calculated in the same manner as for 
our clusters with corrections applied for photomet- 
ric redshift errors and field subtraction. However, 
since the magnitude limit of the SDSS at low red- 
shift probes much deeper down the cluster lumi- 
nosity function, no correction for incompleteness 
need be applied for measuring the luminosity func- 
tion to Mg + 3. The luminosity function for each 
cluster is reported in Table [3] within a selection 
radius of i?2oo- 



A large dispersion in the value of M B and a 
is present for the clusters. Averaging all of the 
clusters together, we find M B = —20.55 ± 0.56 
and a = —0.84 ± 0.32 within -R200, which is 
a significantly larger dispersion than if we had 
first summed the clusters together and then mea- 
sured the luminosity function (see Appendix lAl for 
further issues with ensemble-averaged luminosity 
functions) or from the measurement error associ- 
ated with an individual cluster (typically around 
a a ~ 0.15). For R = 1 Mpc, we find values of 
M* B = -20.39 ±0.48 and a = -0.71 ±0.32. These 
values are very comparable to similar studies once 
converted to our magnitude system. Figure [12] 
and Table 2] contain values and references to other 
measurements of M B and a for low redshift clus- 
ters. 

4. Evolution of the RSLF 

The central question of this is work is deter- 
mining changes in the cluster RSLF shape with 
redshift. To this end, the luminosity function pa- 
rameters, M B and a, are plotted in Figure [12] for 
individual clusters in our low- and intermediate- 
redshift samples, their mean values binned in red- 
shift, and other clusters' values published in the 
literature. The latter have been transformed into 
our magnitude system. Measurements of the lumi- 
nosity function based on isolating the red-sequence 
via morphology, single-function fits to deep lumi- 
nosity functions (M* + 5; see §3.3 and 5.2), and in 
galaxy groups have been excluded. To provide the 
closest comparison to other studies, we plot the 
data for a cluster radius of 1 Mpc. Values for the 
parameterizations within i?2oo instead of 1 Mpc 
generally are closer to a = — 1 with a brighter M* 
for both the low and intermediate redshift clus- 
ters in our sample, but give qualitatively the same 
trends with redshift. 

Our clusters exhibit an increase in M B with 
redshift, as also seen for the field RSLF (Willmcr 
et al. 2006). However, the cluster RSLF M* is 0.5- 
1 mag brighter in the field value at 0.7 < z < 0.9, 
but 0.5 mag fainter at z ~ 0.5-i.e., brighten- 
ing with redshift appears steeper in clusters for 
z > 0.5. In contrast, the value of cluster RSLF 
M* is relatively flat between z = — 0.5, consid- 
ering in concert local values from our low redshift 
sample or the literature. Overall, the brightening 
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of rest-frame -B-band M* in clusters and the field 
agree well with the findings by de Propris et al. 
(1999) and Lin et al. (2006) in the near- infrared 
for the fading of simple stellar populations formed 
between 1.5 < z < 3. 

We find no trend in the faint end slope, a, with 
redshift within our sample. For the low redshift 
clusters, a shows a large range of values. How- 
ever, the average value for the low redshift clus- 
ters within R = 1 Mpc is a = -0.68 ±0.42, which 
increases to —0.71 ± 0.29 if we only include clus- 
ters with a v > 1000 km s . Using this number 
for the high mass clusters (which is probably most 
comparable to our cluster sample), only one clus- 
ter, MS0451, is significantly different at the 1. To- 
level. Including this cluster, the intermediate red- 
shift sample has an average of a = —0.59 ± 0.26, 
which is not significantly different than the aver- 
age value for the low-redshift cluster sample. 

To further investigate the question of evolution, 
we plot a as a function of cluster velocity disper- 
sion for our full sample of low and intermediate 
redshift clusters in Figure 13. The data are cal- 
culated within a cluster radius of IMpc and i?2oo- 
The mean of the two distributions is different but 
not at a significant level. From our sample of inter- 
mediate redshift clusters, we find no significant ev- 
idence for a deficit of galaxies occurring to M* + 3 
as compared to low redshift clusters. 

5. Discussion 

5.1. Assessment of Results in the Litera- 
ture 

Recent literature contains claims and refuta- 
tions of a deficit of faint, red galaxies in intermedi- 
ate redshift clusters. A deficit would be significant 
because a change in the RSLF shape implies an on- 
going or multi-epoch formation scenario, beyond 
passive evolution of a coeval population forming 
at high redshift. We suspect different conclusions 
regarding the deficit (or lack thereof) arise in part 
from differences in samples, analysis methods, or 
comparisons to local samples. We enumerate these 
points below. For comparison and reference pur- 
poses during this discussion, we list all of the rele- 
vant studies of the red sequence cluster luminosity 
function and their principal attributes in Table [4j 



5.1.1. Analysis methods 

A number of studies measure the luminosity 
function shape in terms of a dwarf-to-giant ratio 
(DGR). Typically, the DGR is defined as a ratio 
between the number of galaxies in two magnitude 
bins, and as such, it is a much simpler calcula- 
tion than the luminosity function. However the 
DGR suffers from a number of problems rarely ad- 
dressed in the literature. First and foremost, the 
DGR has had a range of definitions (c/. Ferguson 
& Sandage 1991, Seeker & Harris 1996, Driver, 
Couch, & Phillipps 1998, Tanaka et al. 2005, de 
Lucia et al. 2006, Gilbank et al. 2007, Koyama et 
al. 2007, Stott et al. 2007 ), and consequently is 
not directly comparable between all studies. Often 
the DGR is based around observational (i.e., de- 
tection) limits such that the definitions of "dwarf" 
and "giants" do not necessarily reflect any natural 
division between galaxies and ignores traditional 
splits between such systems. A further weakness 
of the DGR is that its value is dependent on ac- 
curate measurements from each of two magnitude 
bins, whereas the luminosity function (as shown 
in Section [3. 2|) is relatively insensitive to variation 
in a single bin. As the faintest magnitude bin is 
likely to have the largest uncertainties due to in- 
completeness, photometric errors, redshift errors, 
or other sources, an accurate statistical calculation 
of the DGR has far greater errors than a similar 
value derived from the fit of the luminosity func- 
tion. 

One of the purported strengths of the DGR is 
that it does not suffer from the covariance be- 
tween M* and a that is a common complaint of 
Schechter-function fitting. However, without mea- 
suring the shape of the luminosity function, it is 
difficult to ascertain whether changes in the DGR 
are due to changes in the bright end (dominated by 
M*) or changes in the faint end (dominated by a). 
While the DGR is strongly correlated to changes 
in a, it is also weakly correlated with changes in 
M* . We illustrate this in Figure [TT1 Following the 
definition of de Lucia et al. (2006) for the DGR: 

DGR = N(M V < -20.0)/iV(-20.0 < M v < -18.2), 

(2) 

it can be shown the DGR will vary approximately 
as logio(DGR) ~ —0.67a for constant M* and 
logio(DGR) ~ 0.62A/* for constant a. Variation 
in the DGR therefore may not only be due to evo- 
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lution in a. As will be discussed later (also see Ap- 
pendix A for issues with ensemble-averaged lumi- 
nosity functions) cluster-to-cluster variation and 
sample selection may affect the measurement of 
the DGR to bias the results and confuse variation 
in the bright-end normalization for evolution in 
the faint end slope. For these reasons we prefer a 
parametric estimate of the luminosity function via 
the Schechter fitting-function. 

5.1.2. Survey quality 

Significant variations exist between surveys in 
terms of the level of analysis and the quality of 
the data. For example, a significant difference 
between our study and most RSLF studies in 
the literature is the amount of simulations un- 
dertaken to understand the completeness and se- 
lection functions of the data. Few other studies 
have undertaken extensive analysis of their incom- 
pleteness, with several notable exceptions. For 
example, Barkhouse et al. (2007) undertook a 
number of simulations outlined in their appendix 
to understand the effects of projection, Edding- 
ton Bias, and color selection. Mercurio et al. 
(2006) performed extensive completeness simula- 
tions for analysis of a low redshift cluster, but 
then adopted conservative luminosity limits well 
above their detection- and selection-completeness 
limits. Andrcon (2006) presented an in-depth 
description of the statistical method for calcu- 
lating the luminosity function in the limits of a 
background population. However, of these exam- 
ples, the first two focus on low-redshift clusters. 
Only Andreon's (2006) study is at comparable 
intermcdiatc-rcdshifts as our own, and, notably, 
finds similar results as our study. 

Imaging depth also varies substantial from sur- 
vey to survey. None of the intermediate redshift 
studies explore beyond M* + 3.5, and many are 
substantially (1-2 mag) shallower, in contrast to 
many of the low-redshift studies, Shallow surveys 
are not able to probe the faint end of the luminos- 
ity function, and are really just probing the giant 
population. 

Another salient difference between surveys is 
wavelength coverage, and the impact this cover- 
age has on reliable selection of red-sequence galax- 
ies. Our study, although smaller in the number 
of clusters than many of the other studies, has 
far greater wavelength coverage and information 



to constrain and identify the cluster populations. 
Many surveys use only two bands and background- 
subtraction to identify cluster sources. In con- 
trast, it has been shown that multi-band data suf- 
ficient to construct precision photometric redshifts 
can be far more reliable in selecting cluster galax- 
ies (Brunner & Lubin 2000, Rines & Geller 2008). 

5.1.3. Band-pass effects 

As has been shown by Goto et al. (2002), the 
slope of the RSLF is passband-dependent. We 
have avoided comparing our results to any of the 
papers that measure the K-band luminosity func- 
tion, but the conversions between other optical 
bands may still suffer from issues beyond sim- 
ple color transformations. Indeed, further com- 
plications between measuring the RSLF in differ- 
ent passbands was initially shown by Smail et al. 
(1998). When measuring the apparent I-band lu- 
minosity function for galaxies in z ~ 0.25 clusters, 
a selection in apparent (U-B) resulted in a deficit 
of galaxies as compared to the same sample but 
selected in (B-I), which covered the 4000 A break 
at that redshift. To maintain consistency, An- 
dreon (2008) always used filters saddling the 4000 
A break regardless of cluster redshift. However, 
a number of studies (e.g., de Lucia et al. 2006, 
Gilbank et al. 2008) have not done this. 

5.1. 4- Radius effects mixed with sample selection 

A further difficulty in comparing the luminosity 
functions from different studies is the selection ra- 
dius adopted in analysis. The luminosity function 
does show variations with cluster radius (Popesso 
et al. 2006, Barkhouse et al. 2007, and §3.1.1 
here) ; a selection of a small, fixed radius will bias 
results for clusters over a range of masses. This 
affect can be especially pronounced when compar- 
ing sources at different redshifts. The Stott et al. 
(2007) sample is an example of the possible prob- 
lems that may arise when using a small selection 
radius of fixed metric size. The low redshift clus- 
ters in their sample have smaller mass (in their 
case, they use x-ray luminosity as a mass proxy) 
than the higher-redshift sample. However, they 
use a constant measurement radius of 0.6 Mpc for 
all of their clusters. For the intermediate redshift 
clusters, they are only measuring the core of the 
cluster where the decline in the faint-end of the 
luminosity function tends to be more pronounced, 
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but for the low redshift clusters, they measure out 
to -R200 ■ Unsurprisingly, they measure a deficit for 
the higher-redshift clusters. Indeed, for the one 
low redshift cluster of similar x-ray luminosity as 
the intermediate redshift clusters (Abell 3555 in 
their sample with a DGR = 1.89), they find a 
DGR similar to the intermediate redshift clusters. 
In contrast, Andreon (2008) also uses a nearly con- 
stant radius of 0.5 Mpc for the high redshift clus- 
ters in his sample, and he finds no trends with 
redshift, but fairly significant scatter between clus- 
ters. 

5.1.5. Sample differences 

Our selection of clusters samples the most mas- 
sive structures in the Universe - clusters that are 
almost of factor of two more massive than the 
EDisCS sample (de Lucia et al. 2006), as shown in 
Figure Q3] According to the simulations of Wech- 
sler et al. (2002), the typical line of sight velocity 
dispersion of our clusters, today, should be around 
cr = 1500 km s _1 , placing their masses well above 
the Coma cluster. None of the local comparison 
surveys, or even the CNOC cluster sample (Yee 
et al. 1996) at lower redshift are a fair match in 
mass to our sample; the large sample of Popesso 
ct al. (2004, 2005, 2006) with velocity dispersions 
estimated from their x-ray luminosity are also far 
less massive for the most part. The closest com- 
parison is the Smail et al. (1998) sample with a 
few clusters overlapping our sample, but with a 
much lower average mass. The average mass of 
the Stott et al. (2007) intermediate mass sample 
is only slightly less than our sample, but their low 
redshift sample is far less massive. (We estimate 
the velocity dispersion for the Stott et al. sam- 
ple either from their x-ray luminosities, or from 
the velocity-dispersion listed in Andreon (2008) 
for the overlapping subset.) In Andreon (2008), 
the very high redshift clusters are, on average, 
a factor of 30% less massive than the z ~ 0.55 
sample. The masses of clusters in the Barkhouse 
et al. (2007) and Gilbank et al. (2007) studies 
are estimated using the B gc parameter, which is a 
measure of the cluster richness. As they mention, 
this estimator may have dependencies related to 
redshift. According to their conversions, most of 
their sample is far less massive than those stud- 
ied here. Our low-redshift clusters do span a large 
mass range, but a third of the clusters do have 



a v > 1000 km s _1 , which is comparable to the 
predicted masses of our cluster sample. If lumi- 
nosity function shape, or its evolution, depends on 
cluster mass, than cluster sample selection could 
be responsible for the different results found in the 
literature. 

For illustration of possible sample effects, we 
plot in the bottom panel of Figure [12] the amount 
of evolution in a (the faint-end slope) seen by Stott 
et al. (2007), assuming the value for a of our low 
redshift clusters is correct. To make this compari- 
son, we have converted their measure of the DGR 
evolution (a power-law dependence in 1 + z) to a 
function of a (see Figure 11) using the relationship 
found between DGR and a in i j5.1.1l and assuming 
constant M* . As can be seen, this trend is incon- 
sistent with our data, even when normalizing to 
our low-redshift data. Fitting to our own data, we 
find a much shallower trend which is consistent 
with no change. It is conceivable, we suggest, the 
difference in redshift-trends are due to differences 
in the two samples' cluster masses, and systematic 
trends in cluster-mass with redshift in the Stott et 
al. sample. 

In support of this argument, Koyama et al. 
(2007) find a trend between cluster mass and a in 
the sense that more massive clusters have steeper 
turnovers. We see a similar trend in our inter- 
mediate redshift clusters (Figure 13). However, 
the strength of the trend appears sensitive to the 
radial selection, and our statistics are poor. In 
contrast, there is no apparent trend in our low 
redshift data. Andreon (2008) also sees no trend 
using a larger sample of clusters. On balance, 
while cluster-mass may play a role in explaining 
differences between survey results on the RSLF, it 
is more likely that cluster-mass differences within 
samples, particularly when correlated with red- 
shift, plays a more substantial role in driving ap- 
parent evolutionary trends seen in some studies. 

5.2. Challenges to Measuring Evolution 

In the previous section our focus was on iden- 
tifying survey differences that could explain, in 
part, the discrepant results found in the litera- 
ture on RSLF evolution. Here we turn instead 
to the astrophysical nature of the RSLF, and how 
its complexity and variation presents fundamental 
challenges to its measurement. 
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In the low redshift Universe, an over-abundance 
of "true" dwarf, early-type galaxies (those with lu- 
minosities below M B < —17) are found in clusters 
compared to the field (Driver et al. 1994, de Pro- 
pris et al. 2003) and lower density environments 
(Trentham et al. 2005), indicating some relation- 
ship between the environment and shape of the 
faint RSLF must exist. In fact, since the RSLF 
is a composite population (e.g., massive ellipti- 
cals, intermediate mass lenticulars, and low-mass 
spheroidals) a single Schechter function is likely an 
inadequate description of the RSLF; bumps and 
dips have long been noticed in cluster luminosity 
functions (Kashikawa et al. 1995, Jerjen & Tam- 
mann 1996). Since these different sub-populations 
of the red-sequence have different relative densi- 
ties with environment, changes in the shape of the 
RSLF with environment is to be expected at all 
luminosities. This is likely what drives the ob- 
servation that cluster RSLFs changes shape with 
selection radius, and perhaps also differences be- 
tween clusters. 

Due to the composite nature of the RSLF, a 
two-function fit (cither two Schechter functions 
or a Gaussian for the luminous component plus 
Schechter) is commonly used for parameterizing 
low-redshift cluster luminosity-distributions (Bi- 
viano et al. 1995, Popesso et al. 2006). When only 
a single function is fit to the luminosity function, 
the faint end slope is strongly affected by the lim- 
iting magnitude. In intermediate-redshift cluster 
studies, single-function fits are made exclusively, 
yet as Table 4 shows, there is substantial varia- 
tion in sampled depth. 

Since the SDSS data has much greater depth 
then our intermediate-redshift sample, we can ex- 
plore the amplitude of the effect of measuring the 
RSLF to different depths for a given selection ra- 
dius. At M B +5 the faint-end slope should be dom- 
inated by the dwarf population (Trentham & Tully 
2002, Popesso et al. 2006). If we measure a to 
M B + 5 instead of M B + 3, we find a faint end slope 
of a = —1.15 instead of -0.84, i.e., the luminos- 
ity function appears to be slightly rising instead 
of falling at these greater depths. (This value is 
possibly an underestimate as incompleteness will 
start to effect the number counts at Mg + 5 in the 
SDSS data; no correction for incompleteness has 
been applied to these data.) In other words, as 
the sampled depth decreases relative to M* (as we 



would expect and indeed has happened in higher 
redshift surveys), the apparent slope of the RSLF 
increases, i.e., an apparent deficit appears. 

The reason for this effect is simply that the 
bright-end of the RSLF is typically characterized 
in shape as a 'bump,' i.e., having a maximum 
near M* . This bright bump is dominated by giant 
galaxies. The apparent decline at lower luminosi- 
ties is substantially modulated by the relative am- 
plitude of the dwarf-galaxy luminosity function, 
but intermediate redshift surveys do not reach 
down to faint enough magnitudes to well charac- 
terize the luminosity function of this population 
itself. Our data for MS0451 is a good example. 
In this case there is some evidence for an up-turn 
in the RSLF at the faintest magnitudes and a dip 
occurring prior to this upturn about 2-2.5 mag 
fainter than M*, indicative of a composite popu- 
lation. However, the faint-end slope of the lumi- 
nosity function we measure down to M* +3.5 is 
dominated by the excess of bright galaxies. 

Survey-depth alone is not a panacea. The above 
conclusions are modulated by variations seen be- 
tween clusters RSLFs and possible correlations 
with mass. There are indications that the RSLF 
M* is brighter in clusters than in the field (de Pro- 
pris et al. 2003, but Hilton et al. 2005 notes a 
few exceptions), and brighter still in more mas- 
sive clusters (Croton et al. 2005, Hansen et al. 
2005). Similarly, for a sample of 10 clusters, Yagi 
ct al. (2002) find the dip in the luminosity func- 
tion around Mb ~ — 18 becomes much stronger in 
higher-mass systems. 

Unfortunately, there is considerable scatter 
that washes out any clear correlation between the 
shape and normalization of the luminosity func- 
tion and cluster mass for large samples (Biviano 
et al. 1995, Popesso et al. 2003, de Propris et al. 
2003 and 2005, Barkhouse et al. 2007). Within 
our low redshift sample, M B does brighten with 
cluster mass within i?200j but the correlation is not 
statistically significant. Variations in the complex 
luminosity-function shape are seen in the nearby 
clusters of A168 (Yang et al. 2004) and Shapely 
super cluster (Mercurio et al. 2005). In our sur- 
vey, variations are most obvious in comparing two 
of our z <~ 0.55 clusters. Even though both are at 
the same redshift and are of similar mass, MS0451 
has a very pronounced bright end as compared to 
the relatively flat bright end of C10016. 
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Cluster to cluster variation is a significant effect 
and a detriment to measuring evolution from a 
small, diverse sample of clusters. The alternative, 
namely, to studying large samples of clusters to 
average over variations also presents a challenge, 
since each cluster must be observed to adequate 
depth to well characterize the luminosity-function, 
as discussed above. 

Along these lines we close with one further cau- 
tion pertaining to large samples. With signifi- 
cant variation in the luminosity function, the en- 
semble average luminosity distribution may not 
be the same as the average of the individual lu- 
minosity distributions. Simulations we have per- 
formed (see Appendix [AJ indicate systematic er- 
rors do result from ensemble averaging, and are 
particularly large for the DGR index. We con- 
clude that in order to confirm the hypothesis of 
evolution of the red sequence, a far more signifi- 
cant sample of intermediate redshift clusters has 
to be observed to greater depth to constrain the 
true dwarf population. These observations must 
be analyzed in a manner consistent with analysis 
of low-redshift cluster samples, and should not be 
ensemble-averaged. 

5.3. Evolution Scenarios 

If we take the results of de Lucia et al. (2006) 
and Stott et al. (2007) at face value, in concert 
with ours, we conclude the red sequence is al- 
ready in place at intermediate redshifts down to 
faint magnitudes (Mb ~ —18) in the most mas- 
sive clusters, whereas the build-up is still occurring 
in lower-mass clusters. Mass-dependent evolution 
would be in concordance with a general trend 
of down-sizing in hierarchical structure-formation 
scenarios, where the most massive structures form 
earliest. Here, this would have to occur on two 
scales: (i) The most massive clusters would be the 
ones to have their red sequences in place first, (ii) 
For any given cluster, the most massive galaxies 
would evolve to the red sequence first. 

Another more subtle possibility may be occur- 
ring. While the hierarchical structure-formation 
scenario addresses directly the assembly and build- 
up of mass, there is the associated notion that the 
pace of star-formation is also set by this build up. 
Specifically we expect to see, and perhaps do see, a 
down-sizing in the co-moving star-formation rate 
such that it is dominated by lower-mass systems at 



later times. Star-formation, however, is stochas- 
tic on a galaxy scale. Further, small-mass bursts 
in large-mass systems can dramatically change the 
colors, providing a blue, if ephemeral, photometric 
icing. Consequently, while the mass build-up may 
indeed be hierarchical, the growth of the RSLF 
as a function of luminosity may be substantially 
modulated by on-going star-formation. If this is 
the case, the evolution of the shape of the RSLF 
(as distinguished from the normalization) may be 
a rather subtle phenomenon, as we observe within 
our own study. 

By looking at the star-formation histories of 
galaxies currently on the red sequence, we can put 
some additional constraints on the RSLF evolu- 
tion. Studies place the formation epoch of red- 
sequence galaxies, (particularly the massive ellip- 
ticals) at z > 2, with only passive evolution since 
that time (e.g., Holden et al. 2004). However, 
elliptical galaxy stellar populations are typically 
found to be far more uniform than SO populations, 
yet SO's are counted in the RSLF, and are a sig- 
nificant - if not dominant - component at bright 
absolute magnitudes in local samples. If a deficit 
of low-luminosity, passively evolving galaxies ex- 
isted at intermediate redshifts, then we would ex- 
pect to see some correlation between stellar ages 
and magnitude in local cluster red-galaxy popula- 
tions. Perhaps SO's are part of this picture. 

Indeed, some SO galaxies are relatively recent 
additions to the cluster environment. They are 
missing in intermediate clusters relative to local 
clusters (Dressier et al. 1997), and almost 40% 
of those that are found at intermediate redshifts 
have spectra indicative of recent star formation 
(Tran et al. 2007). In some local clusters, SO's 
show far greater spread in their stellar histories 
than elliptical populations (Kuntschner & Davies 
1998), and the color-magnitude relationship is far 
tighter for ellipticals than SO galaxies (Bower et 
al. 1992). 

If an evolving SO population is to be responsible 
for a deficit in the intermediate redshift RSLF at 
faint magnitudes, critical then is determining the 
luminosity range over which the SO's dominate the 
RSLF. It's well known that locally the SO luminos- 
ity function is quite broad, and comparable to E's 
(Binggeli, Sandage & Tamman 1988). Further ev- 
idence of their contribution to the RSLF can be 
found in the scatter in the color magnitude dia- 
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gram. Down to Mb = — 16, which is well below 
the magnitude limit of any of the intermediate red- 
shift cluster surveys, the scatter in the CM relation 
for the red population in the nearby Perseus clus- 
ter is less than 0.07 mag (Conselice et al. 2002). 
The small amplitude of this scatter is typically in- 
terpreted to mean that there has been a very small 
range in star-forming histories in this red popula- 
tion. The lack of trend in the amplitude with lu- 
minosity is also indicative of galaxies having com- 
parable range of formation histories. Greater scat- 
ter is seen at fainter magnitudes (Seeker & Harris 
1997) where the dE population begins to domi- 
nate, but these depths have not been probed at 
intermediate redshift. Similarly, Poggianti et al. 
(2001) find that 50 — 60% of galaxies at all mag- 
nitudes have little to no evidence for any star for- 
mation since z <~ 2. 

These lines of evidence indicate one of three 
possibilities: (1) the notion that SO's are young is 
incorrect; (2) the scatter in the color-magnitude 
diagram is insensitive to the age-variations un- 
der consideration; or (3) the younger population is 
well-distributed in luminosity such that they con- 
tribute to scatter (increase it), but not in a way 
that drives a significant trend with magnitude. We 
suspect the latter scenario is most likely. An ob- 
vious next step is to repeat the RSLF experiment 
with adequate imaging or kinematic data to well- 
distinguish ellipticals from lenticulars in a large 
sample of clusters, probing to depths of at least 
M*+3 and preferably to M*+5. 

6. Conclusion 

In this paper, we have explored the red se- 
quence luminosity function in five intermediate 
redshift galaxy clusters. The luminosity functions 
are measured from deep UBRIz plus narrow band 
imaging from the WIYN telescope. Red sequence 
galaxies are identified from their rest frame colors 
and photometric redshifts within selection radii of 
1 Mpc, 0.25xi?2oo, 0.5xi?2oo, andi? 2 oo- Extensive 
simulations are performed to assure the quality of 
the detection, photometric redshifts, and measure- 
ment of the luminosity function. The quality of 
the data is confirmed through the measurement of 
the field luminosity function from the off-cluster 
sample. 

To provide a low redshift comparison sample, 



we also measured the red sequence luminosity 
function in a set of 59 high mass clusters with 
data from the SDSS. The same process for mea- 
suring the luminosity function for the higher red- 
shift cluster was used here with similar definitions 
for cluster galaxies and the red sequence. For 
both sets of clusters, we find comparable luminos- 
ity functions to those fount in the literature for 
previously studied systems. 

We have two primary conclusions concerning 
the RSLF evolution: 

• Mg evolves in a similar manner as the field 
luminosity function and has faded by about 
0.7 mags over the last 6.5 Gyrs. However, 
little evolution is seen between z = and 
z = 0.5 for the massive clusters. 

• The faint end slope, a, shows no indication 
of evolution between our low and interme- 
diate redshift samples. In addition, we find 
no relationship between the cluster velocity 
dispersion, a v , and a for the high mass clus- 
ters. 

In an extensive comparison to measurements 
from the literature, we have two additional con- 
clusions: 

• Selection effects can be critical to the deter- 
mination of any signal of evolution. Clusters 
do show variations with the luminosity func- 
tion in terms of mass and radius, which can 
lead to erroneous conclusions in terms of evo- 
lution, if not carefully accounted. Although 
small, this survey predominately measures 
the RSLF in massive system which seems 
to display different behaviors than low mass 
systems. 

• Significant cluster to cluster variations ex- 
ists, even at a given mass. The disper- 
sion in cluster luminosity-function parame- 
ters, measured for individual clusters, is typ- 
ically an order of magnitude greater than 
the error estimate on those parameters from 
fitting to an ensemble-averaged luminosity 
function. Significantly more work needs to 
be done to better understand these cluster 
to cluster variations. 
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A clear picture of the evolution of the clus- 
ter RSLF remains elusive, yet such a picture is 
critical for understanding the processes that drive 
the growth and transformations of cluster-galaxies 
over cosmic time. It is conceivable that a suf- 
ficiently delineated map of the RSLF evolution 
with time, cluster-mass, and location within the 
local cluster environment can help confirm or re- 
fute predictions of KCDM models on time-scales 
relevant to the assembly of the most massive, viri- 
alizcd systems in the universe. While the obser- 
vational challenge has yet to be met, the theoreti- 
cal models also fall short on definitive predictions 
due to the complexity of the gas-physics (includ- 
ing star-formation) that channels the transforma- 
tion of galaxies from the blue cloud onto the red 
sequence. For example, clean predictions of how 
and when the different sub-populations along the 
red sequence (e.g., E, SO, dE) are formed are yet 
to be had. As such, the red sequence, as inter- 
preted simply as a mass sequence, will likely con- 
tinue as a critical observational foil thrown up as 
a test of hierarchical models. Tracking the fate 
of the blue population of galaxies that cause the 
Butchler-Oemler effect is a key to understanding 
the physics behind the transformative processes 
in clusters; at least some of these systems are 
likely to be the progenitors of the red-sequence 
galaxies. The complexity of the astrophysics will 
likely stymie a definitive observational picture of 
the transformation process, however headway can 
be made by juxtaposing the blue and red popu- 
lations in the context of environment. In future 
work we will investigate the cluster blue-galaxy 
population in this context. 
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A. Ensemble-average Luminosity Function 

An important aspect of the RSLF brought forward by our work so far is the variation between clusters 
seen even at the same redshift and with similar masses. This phenomenon is fairly well documented for local 
cluster samples. A number of studies measure the luminosity function not for individual clusters but instead 
for the coadded luminosity distribution of an ensemble of clusters. Such studies generally either compute a 
straight average of the luminosity distribution, or a weighted average. The primary method for weighting is 
to normalize each cluster by the number of bright galaxies (e.g., Colless 1989). 

Unfortunately there has been no investigation of whether systematics in the inferred luminosity-function 
parameters are introduced by the cluster-averaging described above, particularly relative to the measurement 
of individual clusters. Since there is a natural need to average over clusters at higher redshift where sources 
are fainter, it is important to understand if this averaging will lead to spurious trends with redshift. 

To explore this possibility, we created a number of simulations using two mock catalogs of 10 and 100 
clusters. We apply errors on the mock cluster-galaxy counts that are representative of those for our low and 
intermediate redshift data. For each of these two catalogs, we generate 100 realizations for each luminosity- 
function distribution we describe below. For each catalog realization, we measured the Schechter-function 
parameters (a, M*, and «&*) and DGR (following de Lucia et al. 2007) for each individual cluster, the average 
of all the clusters, and the weighted average of the clusters following the Colless (1989) prescription. 

For an initial simulation test-set, we assume all clusters have the same luminosity function. Unsurpris- 
ingly, in this sample, all three methods return the same results for the luminosity function and DGR. For 
the remainder of the simulations, we assume a distribution of luminosity functions described by Gaussian 
distribution about M* — —20.0 and a — —0.9. Simulations are generated with distribution-widths in each 
of these parameters of 0.05, 0.10, and 0.25 (with appropriate units). The count-normalization of each lumi- 
nosity function, $* is also assumed to have a Gaussian distribution with widths corresponding to variations 
of 20%, 40%, and 100% about the mean <&* (with the requirement that $* > 0). Simulations were carried 
out where each of the parameters (a, M* , and <&*) was varied individually and also in combination with the 
other parameters. When multiple parameters are varied, the same distribution widths were used for a and 
M* , with the corresponding low-, medium-, or high-percentage width for <i>*. For example, when a = 0.10 
for a and M*, the <i>* distribution width is 40%. One further simulation (lowz) uses variations of 0.25, 0.5, 
and 5% in a, M*, and $*, respectively, which closes matches the values found for the low redshift clusters. 

The results for the simulations are presented in Table [5] for realizations including 10 clusters. The first 
column lists the parameter varied (a, M* , or $*); the second column lists the distribution width for that 
parameter. Subsequent columns list the average and standard deviation for a, M*, and the DGR measured 
from the 100 realizations (a) for individual clusters (columns 3-8), (b) the composite luminosity function 
constructed by averaging the cluster-counts together (columns 9-14), and (c) the same for a composite 
luminosity function constructed through a weighted average of the luminosity functions following Colless 
(1989; columns 15-20 here). For each quantity, the standard deviation is the measured deviation about the 
mean measurement across the 100 simulations. Simulations including 100 clusters yielded identical results 
in the mean, but with smaller standard deviations about the mean for each value. This is due simply to a 
better sampling of the luminosity-function distribution. One exception is for the standard deviation of the 
DGR for individually- measured clusters, which are larger for the 100-cluster simulation (see below). 

For measurements of individual clusters, we are able to recover the original Schechter function values of the 
parent distribution even for large spread in the parameterization. In contrast, the measured DGR does not 
behave in a similar manner as the parameterization. For a single cluster with the nominal luminosity function, 
the DGR would have a value of 4.71. However, as the variation in the luminosity function parameterization 
increases, the average value of the DGR tends to increase as well. The increase in the average value of the 
DGR is due to the Poissian-like distribution in the DGR values due to a Gaussian spread in a. The median 
is a better statistic to measure the DGR from an ensemble of individual measurements than the average. 
For individual measurements of 100 clusters with <r Q= o.25, the median value of the DGR is 4.78, while the 
average is 5.00. 
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The results for the composite luminosity functions indicate that averaged luminosity functions tend to 
yield estimates of increasingly negative a as the variation in a increases, whereas the weighted luminosity 
yields the opposite effect. In both cases, we see a small increase in the measured mean value of M*, with the 
increase comparable to the measured dispersion. Neither case shows any change in parameterization with 
changes in $*. For both cases, the DGR is far better behaved (returning a DGR corresponding to the value 
measured for a in each case) , although the systematic deviation in the mean DGR value is no longer coupled 
with the underlying variation in luminosity function parameters. As expected, the composite luminosity 
function masks the intrinsic cluster to cluster variations. 

For the lowz case, which represents the dispersion measured for our low rcdshift clusters, we find very 
different results from the individual, average, and weighted measurements. The individual cluster measure- 
ments do return the input luminosity function as expected, but with significant variation in a and M* . 
The average DGR is large, but the median, once again, provides a much more accurate statistic. However, 
the average and composite luminosity functions both perform much worse. The average luminosity function 
measures a far steeper slope of a = —1.11 rather than the input value of a = —0.90. The weighted luminosity 
function performs slightly better, but the measured a also decreases to —0.97. Even worse, the measured 
M* increases to M* = —20.48 for both cases (nearly an 0.5 mag shift). Because of this shift in M*, the 
resultant DGR value is generally found to be lower than the nominal value. 

These simulations indicate systematics effects are introduced by coadding data from an ensemble of 
clusters. These effects are minimized by avoiding the DGR formulation and measuring Schechter-function 
parameters. In the case where luminosity functions can be measured for individual clusters, this is clearly 
preferable. When this is not the case, we suspect a maximum likelihood approach that simultaneously fits 
the cluster ensemble assuming a distribution in luminosity function-parameters may be promising. 
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Table 1 

Intermediate- Redshift Cluster Properties 
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DEC 
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rlberg ct al. 19E 


>8, Tran ct al. 


1999, Lubii 


i et 


al. 2004, 


Gal & Lubin 


2004. 





Rcf: Donahue ct al. 2003, Henry 2004, Lubin ct al. 2002, Gioia ct al. 2004 



Table 2 

Intermediate-Redshift Cluster Red Sequence Luminosity Functions 
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Table 3 

Low Redshift Cluster RSLF within i? 2 oo and M* b + 3 
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Table 4 

Summary of Extant Cluster RSLF Studies 
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a For studies with a single cluster, gives the redshift of the cluster. For studies with multiple clusters, gives the range and median value of the rcdshifts. 
b The Field of View (FOV) is in square arcminutcs 

c imiting magnitude in the R band. Transformations are made to the R-band from other optical bands if necessary. Other magnitudes indicated when necessary. All magnitudes transformed 
to Vega. 

^Background refers to the method used for the seperation of field and cluster sources. The diferent techniques are: 'BGS' is background subtraction, 'SZ' is spectroscopic identification, 'PZ' 
is photometric redshift identification, and 'Morph' is using morphological information to identify cluster objects. 

c Combine refers to how the clusters were combined. The different methods are: 'Ind' means the luminosity function was measured for each individual cluster, 'Ave' means a single luminosity 
function fit to the average or sum of all of the clusters, and 'Wei' means a single luminosity function to a weighted average of the clusters. 

* Type of measurement made to assess the density of objects. 'LF' is luminosity function. 'DGR' is for dwarf to galaxy ratio, but the definition differs between the different studies. Before 
each definition, the band (rest or observed) that the measurement is made in is labelled. 

s Tanaka et al. use a density estimate to define clusters vs. field but the estimate corresponds to the viral radius, which is approximately equal to i?200 • Mercurio ct al. also use a density 
estimate to differentiate between low and high-density regions. 

h SDSS magnitudes. 

J HST studies where V, i, and z arc the F606W, F775W, and F850LP bands, respectively. For Stott ct al. 2007, the HST bands for V and i arc F555W and F814W. For Andreon 2008, there 
is only two HST bands that bracket the 4000 Abrcak. Filters used include F555W, F606W, F775W, F814W, and F850LP. 

JMci ct al. 2006, ApJ 639, 81 

^Mei et al. 2006, ApJ 644, 759. Observations of a supercluster located at z=1.265. The two clusters are separated by 2 Mpc and have redshifts of z=1.26,1.27. 



Table 5 

Simulations of Luminosity-Function Measurement Variations for a 10-Cluster Ensemble 
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ct refers to the amount the parameter is being varied in the simulation. 
c "lowz" is a simulation with variations matching our low-redshift cluster sample. The simulation has variations in a, M* , and <i>* of <r a = 0.25. a M* = 0.5, and cr^* = 5%. 




Fig. 1. — R-band completeness profiles as function of object size for the six fields observed for this program. 
Curves (dotted, dashed, solid) are measured for real objects selected from the data with apparent half-light 
sizes of approximately 1.5',1', and 0.5', respectively. For reference, the approximate location of M* for the 
evolving field RSLF (Willmer et al. 2006) is shown as a vertical line. 
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Fig. 2. — Photometric redshift vs. spectroscopic redshift for the full sample of galaxies with spectroscopic 
redshifts in our survey fields. Blue and red galaxies arc defined in terms of rest-frame properties as specified 
in text. 
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Fig. 3. — Histogram of the difference between spectroscopic and photometric redshift as a function of z for 
the sample in Figure 1. The top plot is for all blue galaxies and the bottom plot is for red galaxies, as defined 
in text. 
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Fig. 4. — Photometric redshift error as a function of S/N for measurements in this study. S/N is defined 
within the seeing-matched apertures described in §2. For the sample in Figure 1, the squares represent red 
galaxies; triangles represent blue galaxies. Light gray objects are similar measurements made from the SDSS 
sample using our photometric redshift technique. Photometric redshift error begins to increase dramatically 
below S/N=10, but sources at or above the 50% detection limit have S/N > 20. 
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Fig. 5. — Field RSLF at low (top panel; 0.3 < z < 0.55) and high (bottom panel; 0.7 < z < 0.95) redshift 
in our survey fields. Dashed and dotted vertical lines represent the estimated 50% and 90% detection- 
completeness limits for our data in each redshift bin. Light, open symbols represent raw (uncorrected) 
measurements in our images for redshift regimes free from over-densities (as defined in text). Dark, open 
symbols are corrected for detection completeness. Filled symbols are corrected for both detection com- 
pleteness and cluster contamination due to photometric redshift efforts, as described in text. Unlike cluster 
luminosity-function derivations, here we extend the data beyond where contamination corrections are greater 
than 10% for illustration purposes on the quality of the correction, which is based on simulations. For com- 
parison, the field RSLF measured from the DEEP2 data set (Willmcr ct al. 2006) is plotted as a solid line. A 
reduced-x2 statistic between our data and the DEEP2 Schecter-function (no degrees of freedom) is indicated 
in each panel. 
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Fig. 6. — (left) Rest-frame color-magnitude diagrams for each of the clusters. Only galaxies with photometric 
redshifts within ±0.05 of the cluster redshift and within i?200 are plotted. The redshift selections are 
illustrated in the right panels. Solid lines are the color-magnitude relationship from Willmer ct al. 2006; 
dotted lines are the demarcation between red and blue galaxies (Equation 1 in text), (right) Apparent B — R 
color vs redshift distribution for each cluster field. The solid vertical lines represent the redshift selection 
limits for each cluster, but all sources identified as galaxies in the field are plotted. The expected colors as a 
function of redshift for an non-evolving ellipitical galaxy (solid), Sb galaxy (dotted), and NGC4449 (dashed) 
are overplotted. 



32 




Fig. 7. — (left) RSLF within -R200 m each of our five clusters. Open squares represent raw counts; solid 
squares are the counts after all corrections (see text). The best- fit luminosity function are solid black lines, 
with 68% confidence-limits as gray dotted lines. The best-fit \1 and degrees of freedom is given in each 
frame. The 50% and 90% completeness limits are plotted as heavy dashed and dotted lines, (right) Error 
ellipses (65% and 95%) are plotted for a and M* based on the x 2 distribution. The best-fit value for data 
down to the 50% and 90% detection-completeness limit are plotted as solid and open circles, respectively. 
Best-fit values for ten realizations of the photometry with \[2 larger errors are plotted as plus symbols (see 
text). Gray-scale represents results from jack-knife estimates of the errors. We also plot the best-fit value 
from the low rcdshift cluster sample within i?2oo- 
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Fig. 8. — (left) Red-Sequence cluster luminosity function for each intermediate-redshift cluster as measured 
by different selection radii. Points represent corrected counts within i?200 horn Figure 7; the solid line is 
the best-fit Schechter-function luminosity function to those points. Other lines represent best-fits to data 
extracted with different selection radii: green dot-dashed (0.25i?2oo)> red dashed (0.5i?2oo), blue dotted (1 
Mpc). (right) Error ellipses (68% C.L.) for the best-fit M* and a parameters of the luminosity function 
for all selection radii are shown in the left-hand panels (same line- types and colors). Other measurements 
are shown for C10016 (Tanaka et al. 2005,square) and MS1054 (Goto et al. 2005, hexagon; Andreon 2006, 
triangle). Colors are only available in the electronic version. 
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Fig. 9. — Completeness and reliability of red cluster-galaxy selection due to effects from photometric redshift 
errors. Solid lines represent the fraction of simulated cluster galaxies that have photometric rcdshifts within 
the selection window for the cluster. Simulations are described in text. In the worse case, corrections for 
incompleteness due to photometric-redshift errors are 40% at the 50% detection-completeness limit, and typ- 
ically are only 20-25%. At the 90% detection-completeness limit, the correction for selection incompleteness 
is typically below 15%. Dotted- lines represent the fraction of galaxies in the cluster luminosity function that 
are expected to be contamination from the field. Due to the much larger red galaxy luminosity function 
normalization, only a small fraction (<5%) of field galaxies are expected to contribute even at the faintest 
magnitude bins. 
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Fig. 10. — Apparent size-magnitude locus for the red cluster objects is well offset from the 50% (dashed 
line) and 90% (solid line) completeness limit except for the highest redshift clusters (MS1054, and C11604). 
In these clusters, we may be missing a portion of the dimmest objects which would bias our luminosity 
functions towards a deficit. However, since our completeness corrections include corrections for object size 
and magnitude, this bias is diminished. 
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Fig. 11. — Dwarf-to-giant ratio (DGR) relationship with Schechter luminosity-function parameters a (labels 
at top and x axis) and M* (labeled curves at right). The DGR is primarily a function of a, but is also 
affected by the value of M* , and hence DGR measurement interpretation is degenerate between variations 
in both these quantities. 
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Fig. 12. — Evolution of RSLF M* (top) and slope (bottom) with redshift. Individual data from our survey 
is marked by filled, gray circles and diamonds, whereas the average values are marked by dark circles and 
diamonds (see key). Other studies are plotted as given in the key. All except two field studies (Madgwick 
et al. 2002, Willmer et al. 2006) are for cluster samples. Dashed curves in top panel show evolution in 
Mb of a solar- metalicity, Salpter-IMF simple stellar population with formation epochs of z = 1.5,2,3 (top 
to bottom), referenced to Madgwick et al. at z = 0. The dashed curve in the bottom panel represents 
the observed evolution in the DGR relationship found by Stott et al. (2007) normalized to our low-redshift 
measurement of a (see text). M* rises much more smoothly with redshift in the field (Madgwick et al. 2002, 
Willmer et al. 2006) than in the cluster sample, whereas a shows little to trend with redshift. 
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Fig. 13. — Faint end slope, a, as a function of cluster velocity dispersion for low and intermediate redshift 
clusters. The low-redshift data is from measurements based on SDSS spectroscopic measurements of the 
Abcll sample for magnitudes limits of M* + 3. In the left panel all data are measured within a selection 
radius of 1 Mpc. The right panel shows the same, measured within a selection radius of R,2oo- 
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Fig. 14. — Cluster velocity dispersions (masses for constant size) as a function of redshift. Our smaller 
sample of clusters are far more massive than either the CNOC sample (Yee et al. 1996, Carlberg et al. 1997; 
open shields), the EDisCS samples (White et al. 2005, de Lucia et al. 2006; open diamonds), or even Coma 
(Colless & Dunn 1996; hlled square). Other points are the same as in Figure [T2l Predictions for the growth 
of halos from Wechsler et al. (2002) are over plotted as dashed line for halos with a — 1500, 1000, 500 km 
s _1 at z = 0. 
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